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ABSTRACT 

The discovery of kilohertz quasi-periodic oscillations (kHz QPOs) in low-mass X-ray binaries 
(LMXBs) with the Rossi X-ray Timing Explorer has stimulated extensive studies of these 
sources. Recently, Osherovich & Titarchuk suggest a new model for kHz QPOs and the related 
correlations between kHz QPOs and low frequency features in LMXBs. Here we use their results 
to study the mass - radius relation for the atoll source 4U 1728-34. We find that, if this model 
is correct, 4U 1728-34 is possibly a strange star rather a neutron star. 
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Introduction 



With the Rossi X-ray Timing Explorer kHz QPOs have been discovered in about 20 neutron star 
LMXBs (see van der Klis 1998 for a review). In many cases, two simultaneous kHz peaks (~ 200 — 1200 
Hz) are observed in the power spectra of the X-ray count rate variations, with the separation frequency 
Av roughly constant (e.g., Strohmayer et al. 1996; Ford et al. 1997; Wijnands et al. 1997). Sometimes 



a third kHz peak is detected in a few atoll sources during type I X-ray bursts at a frequency fb equal to 



the separation frequency Ai^ of the two peaks or twice that (see Strohmayer, Swank, & Zhang 1998 for a 
review). These two results suggest that a beat-frequency mechanism is at work, with the upper kHz peak 
at the Keplerian orbital frequency at the inner edge of the accretion disk around the neutron star, the third 
peak at the neutron star spin frequency (or twice that), and the lower kHz peak at the difference frequency 



between them (Strohmayer et al. 1996; Miller, Lamb, & Psaltis 1998) 



Within the beat-frequency model, the burst oscillations Vh and the kHz peak separations Av are both 
believed to be close to the neutron star spin frequency and thus Aiy should remain constant. However, 



recent observations have challenged this interpretation: in Sco X-1 (van der Klis et al. 1997), 4U 1608-52 



(Mendez et al. 1998), and 4U 1735-44 (Ford et al. 1998), Ai/ slowly decreases as the frequencies of both 



QPOs increase. Moreover, in 4U 1636-53 the frequeny t^b of the burst oscillations (or half this value) does 
not match the frequency difference Av between the kHz QPOs (Mendez, van der Klis, & van Paradijs 
1998| ). In the atoll source 4U 1728-34, Mendez & van der Klis ( |1999D find that Av is always significantly 
smaller than z^b, even at the lowest inferred mass accretion rate, when Av seems to reach its maximum 
value, and that Aiy decreases significantly as the frequency of the lower kHz QPO increases. 
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A different approacli was recently invoiced by Osherovich & Titarcliuk; (1999a) and Titarcliuk 



Oslierovich ( 1999 ), wlio proposed a unified classification of kHz QPOs and the related low frequency 
phenomena. In this model, kHz QPOs are modeled as Keplerian oscillations under the influence of the 
Coriolis force in a rotating frame of reference (magnetosphere) . The frequencies of the upper kHz QPO 
branch hold a hybrid frequency relation with the Keplerian frequencies t'k (referred to the lower kHz QPO 
branch): = + (2t'i„)^j where is the rotational frequency of the star's magnetosphere. For Sco 
X-1, the QPOs with frequencies ~ 45 and 90 Hz are interpreted as the 1st and 2nd harmonics of the lower 
branch of the Keplerian oscillations (Osherovich & Titarchuk 1999a ). The same interpretation is applied 
to the 35 Hz QPOs in the atoU source 4U 1702-42 (Osherovich & Titarchuk [l999b|). The observed low 



Lorentzian frequency in 4U 1728-34 is suggested to be associated with radial oscillations in the boundary 
layer of the disk, whereas the observed break frequency is determined by the characteristic diffusion time of 
the inward motion of the matter in the accretion flow (Titarchuk & Osherovich 1999). Predictions of this 
model regarding relations between the QPO frequencies mentioned above compare favorably with recent 
observations of Sco X-1, 4U 1608-52, 4U 1702-429, and 4U 1728-34. 

The discovery of kHz QPO features allows to probe not only the interaction between the accretion disk 
and the stellar magnetosphere, but also the structure of the compact stars involved. In this Letter^ based 



on the work of Osherovich & Titarchuk (1999a) and Titarchuk & Osherovich (1999), we study the mass 
radius (Af — i?) relation of the atoll X-ray source 4U 1728-34, and show that it is possibly a strange star. 



2. The mass - radius relation for 4U 1728-34 

Titarchuk & Osherovich ( |1999| ) suggest a model for the radial oscillations and diffusion in the viscous 
boundary layer of the accretion disk in LMXBs. Their dimensional analysis has identifled the corresponding 
frequencies which are consistent with the low Lorentzian and break frequencies for 4U 1728-34 observed 



by Ford & van der Klis ( |1998D , the predicted values for the break frequency related to the diffusion in 
the boundary layer are also consistent with the observed ones for the same source. The presence of the 
break frequency and the correlated Lorentzian frequency suggests the introduction of a new scale in the 
phenomenon. One attractive feature of the model is the introduction of such a scale in the model through 
the Reynolds number. The best flt for the observed data was obtained when 

Ok = (Af/M0)(i?o/3i?s)^/^(j^/364Hz) = 1.03, (1) 

where M is the stellar mass, i?o is the inner edge of the accretion disk, = 2GM/c^ is the Schwarzschild 
radius, and v is the spin frequency of the star. Given the 364 Hz spin frequency of 4U 1728-34 (Strohmayer 



et al. 1996), from Eq. (1) we derive the inner disk radius 

i?o 9(ak/1.03)2/3mi/^ km, (2) 

where m — M/Mq. Since the innermost radius of the disk must be larger than the radius of the star itself, 
this leads to a mass-dependent upper bound on the stellar radius, 

R < 9(ak/1.03)2/3mi/3 km, (3) 

which is plotted in dot-dashed curve in Fig. § The dotted curve in Fig. |l| represents the Schwarzschild 
radius, which presents the lower limit of R. 

Figure |^ compares Eq. (3) with the theoretical M — R relations (solid curves) for nonrotating neutron 
stars given by six recent realistic models for the equation of state (EOS) of dense matter. In models 
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UU (Wiringa, Fiks, & Fabrocini |1988D , BBBl and BBB2 (Baldo, Bombaci, & Burgio |1997D , the neutron 
star core is assumed to be composed by an uncharged mixture of neutrons, protons, electrons and muons 
in equihbrium with respect to the weak interaction (/3-stable nuclear matter). These models are based 
on microscopic calculations of asymmetric nuclear matter by use of realistic nuclear forces which fit 
experimental nucleon-nucleon scattering data, and deuteron properties. In model Hyp (Prakash et al. 



1997), hyperons are considered in addition to nucleons as hadronic constituents of the neutron star core. 



The curve labeled BPAL12 is a very soft EOS for /3-stable nuclear matter (Prakash et al. 1997), giving a 
lower limiting mass and a smaller radius with respect to a stiff EOS. Finally, we consider the possibility that 
neutron stars may possess a core with a Bose - Einstein condensate of negative kaons (e.g. Glendenning 



& Schaffner-Bielich 1999). The main physical effect of the onset of K~ condensation is a softening of the 
EOS with a consequent lowering of the neutron star maximum mass and possibly of the radius. This is 
shown with the curve labeled K~. It is clearly seen that only the upper part of EOSs UU, BBBl, BBB2, 
and BPAL12 are compatible with Eq. (3). 

Another constraint on the mass and radius of 4U 1728-34 results from the requirement that the inner 
radius Rq of the disk must be larger than the last stable circular orbit i?nis around the star, which is given 



by (Bardeen, Press, & Teukolsky 1972) 



Rn 



i?h(3 + Z2- [(3 - Zi)(3 + Zi + 2Z2)]i/2), 



where 



Zi = 1 + [1 - (a/i?h)2]i/3[(i + a/i?h)i/3 + (1 _ a/i?j^)i/3] 
Z2 = [3(a/i?h)2 + Z2]i/2, 



(4) 

(5) 
(6) 



with i?h — GM/c^ and a = lil/Mc (where / is the moment of inertia and fl = 27rz^ is the angular velocity 
of the star, respectively). For a given mass M, from Eq. (2) we can find Rq. From the inequality Rq ^ R,^s 
and Eqs. (4-6), one can calculate the lower limit of a, and from there the lower limit of the moment of 
inertia /. To account for the inflation of the radius caused by rotation, we adopt the formula 



/ = 0.21Mi?2(i _ R^/R)- 



(7) 



proposed by Burderi et al. (1999), to obtain the lower limit on the stellar radius R^ The results are shown 
in the dashed curve in Fig. l|. Combining the constraints R ^ Rq and Rq ^ Rms reveals that the allowed 
distribution of the mass and radius of 4U 1728-34 must lie between the dot-dashed curve and the dashed 
curve. None of the neutron star EOSs extends to this region. Thus it seems to rule out the possibility that 
4U 1728-34 is a neutron star! Including rotational effects will shift the theoretical neutron star M — R 
curves to up-right in Fig. |l|, making the contrast between the theoretical neutron star models and the 
derived mass and radius range of 4U 1728-34 even worse. 

The difficulty of the neutron star model for 4U 1728-34 suggests that it may belong to another type of 
compact objects, that is, strange stars, which are entirely made of deconfined u,d,s quark matter [strange 



matter). The possible existence of strange stars is a direct consequence of the conjecture (Bodmer 1971 



Witten 1984) that strange matter may be the absolute ground state of strongly interacting matter. Detailed 
studies have shown that the existence of strange matter is allowable within uncertainties inherent in a 
strong interaction calculation (Farhi & Jaffe 1984); thus strange stars may exist in the universe. 



^There are actually two solutions of Eq. (7) with given / and M, one serves as the upper limit of R and the other the lower 
limit of R. But the former solution gives an unphysically small radius, and is rejected. 
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Most of the previous calculations (e.g. Alcock, Farhi, & Olinto 1986; Haensel, Zdunik, & SchaefFer 



1996| ) of strange star properties used an EOS for strange matter based on the phenomenological nucleonic 



bag model, in which the basic features of quantum chromodynamics, such as quark confinement and 
asymptotic freedom are postulated from the beginning, though the deconfinement of quarks at high density 



is not obvious in the bag model. Recently, Dey et al. (1998) derived an EOS for strange matter, which has 
asymptotic freedom built in, shows confinement at zero baryon density, deconfinement at high density, and 
gives a stable configuration for chargeless, /3-stable strange matter. In this model the quark interaction is 
described by an interquark vector potential originating from gluon exchange, and by a density dependent 
scalar potential which restores the chiral symmetry at high density. Using the same model (but different 



values of the parameters with respect to those employed in Dey et al. 1998 ) we calculated the M — R 
relations, which are also shown in solid curves labeled ssl and ss2 in Fig. |l|, corresponding to strange stars 
with maximum masses of 1.44 A/q and 1.32 Mq and radii of 7.07 km and 6.53 km, respectively. It is seen 
that the region confined by the dashed and dot-dashed curves for 4U 1728-34 in Fig. |l| is in remarkable 
accord with the strange star models with masses ^ IAMq. Figure |l| clearly demonstrates that a strange 
star model is more compatible with 4U 1728-34 than a neutron star one. This seems to be important 
since the model of Dey et al. ( |l998[) employs quark masses which decrease with increasing density, a 
manifestation of chiral symmetry restoration (CSR), and the parametrization of CSR can be deduced from 
the mass and radius of the strange star fitted. 

The lower limit on the moment of inertia deduced above does not depend on the EOS. In particular 
Eq. (7) is an extrapolation from neutron stars and it may be doubted whether it is valid for strange stars. 
We have checked that the standard calculated moment of inertia of the strange star of mass M ~ 1.1 Mq 
agrees with Eq. (7) with an error of about 20 %, and the numerical factor 0.21 can change to 0.28 (ssl) or 
0.26 (ss2) depending on the model parameters. The moment of inertia calculation of the strange star with 
variable density layers may not in fact satisfy the assumptions of the standard calculation. It is satisfying 
to find that an EOS-independent estimate can already reproduce the inertia to a reasonable accuracy. 



3. Discussion and conclusion 

We have shown that strange stars are more consistent with the properties of 4U 1728-34 compared 



to neutron stars, i/ the model of Osherovich & Titarchuk (1999a) and Titarchuk & Osherovich (1999) 
correctly interprets the high- and low-frequency QPO phenomena in 4U 1728-34. Since our analysis does 
not involve any unknown physical factors such as the magnetic field strength and topology, radiation 
processes, and the mass accretion rate, it seems to provide the first definite evidence for the existence of 
strange stars. It remains to see whether other compact stars in LMXBs would be identified as strange stars 
from observations of their QPO phenomena. If so, there will be very deep consequences for both physics of 
strong interactions and astrophysics. 

Of course our conclusion is subject to the uncertainty in the parameter Ok, which is determined by 
the best fitting of the observed data. Increase or decrease in will enlarge or reduce the allowed region 
of the mass and radius of 4U 1728-34 in Fig. ^ For example, if Ok — 2, nearly all neutron star EOSs are 
compatible with 4U 1728-34. However, according to Titarchuk & Osherovich ( 1999| ), the dependence of 



on Ok is quite strong: = 38024 - 73076ak + 35732ak^ ^ 35732[(ak - 1.03)^ + 0.019]. A ~ 15% variation 
of Ok from the preferred value 1.03 will double the value of x^, greatly reducing the confidence level of the 
fit. Thus we expect that the possible change in ak and in the derived M — R relation is rather small. 
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If 4U 1728-34 is a strange star, Fig. |l| reveals that it possesses a raass I.IM©, which is quite 
small compared to the ^ 1.4M0 masses of standard neutron stars like PSR 1913+16, though still barely 
consistent with the estimated masses of quite a few radio pulsars within 68% confidence limit (see Thorsett 



& Chakrabarty 1999 ). Since the mass accretion rates in LMXB atoll sources are lower than the Eddington 
limit accretion rate by 1 — 2 orders of magnitude, it is conventionally thought that in such systems almost 
all of the matter transferred from the companion stars should be accreted by the compact stars. The light 



mass of 4U 1728-34, however, seems to support the conclusion of Thorsett & Chakrabarty (1999) that there 
is no evidence for extensive mass accretion occurred in LMXBs, implying that even the simplest evolution 
of LMXBs is far from being fully understood. However, formation of strange stars may be different from 
that of neutron stars (i.e. type Il/Ib supernova explosions) concerning the glitching phenomena observed 



in some radio pulsars (e.g. Alpar |1987| ), suggesting that some exotic formation channels may be at work. 
Accreting massive neutron stars in binary systems are not favorable progenitor candidates of strange stars, 
at least in the case of 4U 1728-34, since it requires too large mass ( ^ 0.7Mq) ejected during the conversion 
as the neutron star masses exceed ^ I.8M0. Accretion-induced collapse of white dwarfs remains an open 
possibility of strange star formation, and could be responsible for the light mass of 4U 1728-34. The detailed 
investigation of this subject is beyond the scope of this paper and wil be discussed elsewhere. 

We are grateful to Lev Titarchuk for helpful and prompt e-mails. X. L. was supported by National 
Natural Science Foundation of China. S. R., M. D. and J. D. were supported in part by DST grant no. 
SP/S2/K18/96, Govt, of India. 
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Fig. 1. — Comparison of the M — R relation of 4U 1728-34 determined from RXTE observations with 
theoretical models of neutron stars and of strange stars. The range of mass and radius of 4U 1728-34 is 
allowed in the region outlined by the dashed and dot-dashed curves. The solid curves labeled UU, BBBl, 
BBB2, BPAL12, Hyp, and K~ represent various M — R relations for realistic EOSs of nonrotating neutron 
stars; the solid curves labeled ssl and ss2 are for strange stars. 



